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Abstract 


This report describes the experimental conditions under which the NASA-GSFC 
High Energy Cosmic Ray Experiment was calibrated with the G-10 + 4.7° test beam 
at the AGS at Brookhaven National Laboratory, A detailed discussion is given 
of the momentum of the test beam used. In addition, summaries of the data runs 
taken are included. 

I . The Problem 

Cosmic rays are presently the subject of intensive study because of their 
relevance to the field of astrophysics. Astrophysics is concerned with the nature 
of extraterrestrial phenomena, especially with the nucleo-synthesis processes 
involved in building up the nuclear isotopes. Since primary cosmic rays bring 
us a sample of material from outside the solar system, a study of their charge, 
energy, and direction of arrival yields information on their formation and 
acceleration and on the properties of the interstellar medium through which 
they propagate. 

15 

Cosmic ray energies above 10 eV have been studied mainly by means of 
extensive air showers, and energies below 10^^ eV have been studied extensively 
using balloons and satellites. These studies have produced models of cosmic 
ray acceleration, propagation, and storage which can best be tested by studying 
the energy spectra and charge composition in the energy range 10^*^ to 10^^ eV. 

At these energies the number of surviving primaries is too small relative to the 
secondary flux at sea level, or even at mountain altitudes. Therefore, experi- 
ments are being carried out at high altitudes (125,000 ft.) with equipment 
carried by research balloons, a method that allows the study of energies up 
to lO eV. Present experiments in the Soviet Union and planned experiments 
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in the United States use satellite-borne instruments to extend this energy range 
up to 10^^ eV. 

II . The Instrument 

An instrument to carry out the measurements described above has been designed 

and built at Goddard Space Flight Center by a group including J. F. Ormes, F. B, 

McDonald, and V. K. Balasubrahmanyan of the High Energy Astrophysics Branch. 

This instrument contains four digitized spark chambers to define the 

trajectory of the incoming particle (see Figure 1). It also contains a charge- 

determining module consisting of two plastic scintillators that produce an 

2 2 

output proportional to Z /3 , a CsI(Tl) scintillator that produces an output 
2 2 

proportional to Z /p , and a Cerenkov detector with an output proportional to 
2 1 

Z (1 ~ — 7 — j™) . Together, these four counters enable the determination of the 
P n 

charge Z of the incoming particle up to Z = 26, and they also provide 
triggers for the spark chamber system. 

For energy measurement, an ionization spectrometer (IS) is used. This 
consists of alternate layers of high-Z absorber and plastic scintillator. Each 
scintillator section is viewed by a photo-multiplier tube which registers a 
pulse proportional to the number of particles traversing that section of the 
spectrometer. The cascade shower induced by the incoming particle results in 
many secondary particles. Measurement of the number of secondary particles at 
many points in the shower alloxvs reconstruction of the energy of the incident 
particle. The top absorber layers are made of tungsten to study electrons by 
means of electromagnetic cascade showers, while the lower layers are made of 
iron to study incoming nucleons and nuclei by means of their nuclear cascade 
showers. The electron section contains 12 radiation lengths of tungsten, which 
is only 0.4-4 nuclear mean free paths. This is followed by up to 5,5 mean free 
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paths of iron in the nuclear cascade section. 

The equipment uses one of four triggering modes, depending on the type of 
particle to be observed: electron mode, proton mode, heavy nuclei mode, and 
calibrate mode. The calibrate mode triggers on any particle within the aperture 
of the equipment, while the others are self-explanatory. 

Physically, the ionization spectrometer weighs 3 tons. For flight purposes 

it is contained in a housing 5 feet in diameter by 8 feet long. The geometry 

factor for particles that traverse at least 3 mean free paths of material is 
2 

about 600 cm sr. The electronics readout time (dead time) is 0.05 sec, but 
the detector remains active only for 2 microseconds after passage of the first 
particle. The internal layout is shown in Figure 1. 

III. Calibration Requirements 

In order to use the IS to measure particle energies, it is necessary to 
make a theoretical analysis of detector response as a function of energy. This 
will also be a function of arrival direction and particle type, and is best 
carried out using Monte Carlo techniques. Such calculations have been carried 
by V. Jones at Louisiana State University, and have influenced the design of 
the IS. 

However, in order to make absolute energy measurements, the IS response 
must be calibrated using particles of known energy. This is necessary in order 
to determine the energy going into invisible processes such as nuclear excitation 
as well as the efficiency of various components in the system. Also, since the 
energy of each incoming particle is only sampled by the detector, the accuracy and 
resolution in these measurements depend on the fluctuations in the energy con- 
tained and energy lost by the IS. The distribution of these fluctuations must 


also be studied in the calibration run. 
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Since the instrument is designed for energy measurements in the range 10 

14 

to 10 eV, this calibration has been performed at the highest available machine 
energy. Lower energies were also used in order to allow a good check of the 
theory which is to be used to extrapolate to even higher energies. Because the 
fraction of the incident energy which escapes from the detector depends on the 
angle of incidences the calibration has been carried out at several different 
angles. Calibration of the electron section has previously been accomplished 
using 5.4 to 18 GeV electrons at SLAG. Calibration of the nuclear section is 
the subject of this report. The data for this purpose were obtained using a 
proton beam at the AGS, 

IV. Physical Layout 

12 3 

The experiment was set up in the G-10 + 4,7° Test Beam ’ ’ in the East 
Experimental Area at the AGS. The location was 2258 inches from D 2 S approxi- 
mately the same as the MIT-Brown Experiment (#416). The configuration of the 
beam transport system is given in Table I and illustrated in Figure 2, Because 
of power supply limitations, the maximum available momentum was about 20 GeV/c. 
Of the quadrupoles, only and Qg were available during these runs. Power 
limitations prevented the use of Q 2 , and the power supplies for the other quad- 
rupoles had been removed. A brief description of the experimental apparatus in 
the beam upstream of our experiment is presented in Table II. 

The IS xi?as supported by jackscrews on a wheeled cart, allowing the hori- 
zontal position and orientation as well as the vertical position to be changed 
with respect to the incoming beam (see Figure 3). 

The whole experiment was encased in a large light-tight tent to prevent 
light leakage into the photomultiplier tubes. 
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V. Data Acquired 

The experiment arrived at BNL on 5 July, 1970 and was set up in a vertical 
position for test runs on cosmic ray muons. It was then moved into position in 
the beam and test runs conducted to set up triggering modes and the beam tele- 
scope during the period 13 - 17 July. The production runs were made from 
18 July through 21 July. Following the production runs, a cosmic ray muon run 
was taken with the AGS off and the IS in the beam position. The experiment was 

removed from the beam and left BNL on 25 July, 1970. 

Twenty four production runs of 5 to 10 thousand events each were made, as 

described in Table III. These runs will allow the study of the dependence of 

detector response on; 

(1) Momentum ( from 10 to 20 GeV/c) 

(2) Angle of Incidence (from 0° to 20°) 

(3) Position of Incident Particle (various locations over the aperture, 
see Figure 4) 

(4) Incident Particle (protons and pions) 

(5) Depth of Spectrometer (3.5 to 5.5 nuclear mean free paths). 

These measurements will allow the objectives listed in Section III above 
to be reached, with the exception that the highest momentum used was only 20 GeV/c 
rather than the desired 30 GeV/c. 

VI. Triggering Requirements 

Aside from triggering requirements which are internal to the IS itself, the 
following external requirements were established: 

B Sj^ S 2 S-j N C for protons 

B S^ S^ N C for pions 



- 6 - 


B is a gate which is only open during the flat top of the beam spill. S2 

are scintillation counters in coincidence, forming a beam telescope as described 
in Table I. N is an anti-coincidence trigger fired by noise. C is a signal from 
the gas Cerenkov counter. 

The anti-noise circuit was required to suppress noise pulses generated by 
the many spark chambers on the floor. A wire antenna a few feet long was used 
to pick up the noise, which was amplified (XIO) , sent through a discriminator, 
and finally to the anticoincidence gate. 

The gas Cerenkov counter was filled to twice the threshold pressure for 
kaons at each beam momentum, a value which is about half the threshold pressure 
for protons. This allowed discrimination between protons and other particles 
such that at least 987. of all particles accepted when the Cerenkov is in anti- 
coincidence are protons. 

Under these same conditions, but with the Cerenkov in coincidence, the 
+ 10 GeV/c pion beam contains less than 97. kaon contamination. The + 20 GeV/c 
pion beam, however, requires a different setting of the gas pressure to get 
adequate discrimination against kaons, a procedure which was not followed in 
this particular experiment. The threshold pressures for various particles and 
momenta are listed in Table IV. The Cerenkov counter is built to operate at 
pressures up to 70 psia and above this pressure the safety valve will release. 

VII, Beam Momentum 

One of the most difficult aspects in the analysis of this experiment is 
the determination of the central value of the beam momentum. Because the first 
collimator and the effective target position are not on the magnet axis, as 
shown in Figure 2, and because most of the focussing quadrupoles are not avail- 
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able for our range of beam momentum, the behaviour of the beam is far from ideal, 
A. Momentum Resolution 

In a beam transport system, momentum definition is accomplished by constrain- 
ing the beam to lie along a line, to transit a particular bend angle, and to pass 
through a point. The momentum resolution is then determined by the angular 
resolution due to the finite width of the beam "line" and the aperture "point". 

In this test beam the momentum resolution is determined by the angular 
resolution in the first bend angle as follows: 

(1) The target T^ and the exit aperture of the first collimator 
form a line with 


6 G = 


1.0 


3.50 mrad 


268.25 + 18.0 

(2) The exit aperture of the second collimator C 2 forms a point with 


6 9 = 


0.5 


= 1.03 mrad 


460.75 + 24.0 

(3) Together this gives 6 0^^ = 4.53 mrad, which is the full width measured 
to the limiting rays. This is to be compared to the nominal bend angle 
= 59.3 mrad, thus 


6 0^ 4,53 


59.3 


c 

'1 


= 7.6% 


(4) If one assumes a transmission function of triangular shape (see below), 
then the full width at half maximum is half the width to the limiting 
rays, i.e., 4 0 = 0.5 6 0 and 
4 P _ A 0T 


0 - 


3.8% FWHM 


A 


This expected value for the resolution agrees closely with a value of 


4 % 


obtained directly using a lead glass detector at 8 GeV/c, 
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If the variable collimator is closed to less than 0.58 inch width the 
resolution is somewhat improved. However, there was no monitor available for 
the width of during the experiment so the actual value cannot be calculated. 
Nevertheless it can be seen that at best, closing completely would give 


A P 
P 


= 3% (FWHM), 

As for the second bend, a similar analysis gives 


A 00 

£ = 207o, which is 

02 


obviously not the limiting aperture in this particular configuration. If 
better momentum resolution is desired, one can select only events which lie 
in a small area of the spark chamber in the NASA experiment. This is equivalent 
to using a narrower beam telescope after the second bending magnet, and could 
allow the resolution of the second bending system to become the limiting .factor . 
B. Tuning Procedure 

The tuning procedure followed was similar to that used by previous groups; 

(1) Set D 2 at the nominal value obtained by linear scaling from the 8 GeV/c 
MIT- Brown setup. (1300 amps corresponds to 8 GeV/c). 

(2) Tune and maximum flux in the beam telescope S 2 S 3 . 

(3) Close the variable collimator C 3 until about 100 particles/AGS pulse 
are obtained in $2 S^. Repeat (2) and (3) as necessary. 

The difficulty arises when this method of tuning is followed because the 
momentum calculated directly from the current^ in differs significantly from 
that calculated directly from the current in D 2 , as shown in Table Va. Further- 
more, the current required in is a strong function of the current in , 

confirming that the beam is off-akis and considerable steering is being done by 
the quadrupoles. 
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C . Central Value of Beam Momentum 

The momentum can be calculated from the magnet currents and the bend angle 
in either or D 2 . However, the bend angle in the system associated with 

is uncertain because the beam coming from the target is off-axis (see Figure 2) 
and because an undetermined amount of steering is done by the quadriipoles Q 5 and 
Qg . The bend angle in the D 2 system is also uncertain because of the large 
angular acceptance of the telescope S 2 S 3 . The problem is to determine the 
beam momentum as accurately as possible in the face of these uncertainties. 

Although the data were not recorded for this purpose, a further check on 
the momentum is provided by the proton thresholds on pressure curves taken with 
the gas Cerenkov counter. The momentum values obtained from the Cerenkov data 
are about 10% higher than the corrected values from the magnets, except at 10 GeV/c, 
However, the pressure curves are very steep near the threshold and they were meas- 
ured in rather coarse steps, so the uncertainty is quite large in these deter- 
minations. Therefore, we rely on the magnets for momentum determination, and 
include the Cerenkov data merely for completeness. It would be possible in 
future experiments with this beam to use the Cerenkov counter specifically 
to determine the momentum, provided that the proton thresholds were measured 
within 1 psi or better. In the only instance where our measurements were of 
reasonable precision, at 10 GeV/c, the momentum from the Cerenkov agrees with 
that from the magnets as seen in Table Va. 

D. Momentum from the D 2 System 

In order to approach the evaluation of momentum from the ^2 magnet system, 
it is necessary first to consider the transmission properties of a beam telescope 
with apertures of unequal width. This is discussed in Appendix A, This background 
makes it possible to discuss the beam momentum in terms of the tuning procedures 


actually used. 
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TgCp) 

o 


o 


m 


(1) It is known that the number of protons F(p) 
leaving the target along the beam line decreases 
rapidly with increasing proton momentum. In fact, 

p~^ or steeper in the region in question, 
„-4 




dQdp 
so that F(p) 




(2) It has been shown that the transmission func- 
tion of the second magnet system has approximately 
trapezoidal shape centered about p 2 , with Pj^ and 
p^ at the middle of the rounded shoulders. P 2 
is the nominal momentum for which the current in 


Pl Pe Ph 


D 2 is set. 



(3) When the production spectrum F(p) and the 
second bend transmission function T 2 (p) are folded 
together the result is a curve which is similar 
to a trapezoid with a tilted top. 

(4) The transmission function of the first magnet 
system has a narrow triangular shape. In tuning, 
the current in is adjusted for maximum flux, 
while D 2 is held fixed. That is, the integral 

Jo F(p) T 2 (p)Tj^(p)dp is maximized. A rough calcula 
tion indicates that the maximum occurs when p^^ = 
MOMENTUM , p that is T^^ is centered at the lower shoulder of T 2 . 
The uncertainty in this tuning is taken as + 69c/2’ the width of the rounded 



shoulder. For the particular system described here we then find 


Pi 


/A 


P2 


V P2 


= 10.1 + 2 . 2 % 


- 11 - 


Thus the actual beam momentum will be approximately 10.1% below the value corres- 
ponding to the current set in magnet 

E. Momentum from the System 

Another approach to the determination of beam momentum is to consider the 

situation at the first bending magnet. The bend angle at is nominally 59.3 

mrad, but because of the fringe fields of the G-11 ring magnet the effective 

target point for positive beams lies outside the beam line. Also, the first 

collimator is displaced approximately 0.64 inches parallel to the beam line 

so that it defines an actual beam line somewhat different from the nominal one. 

The offset of the virtual target Ty varies from l.CDin at 20 GeV/c to 1.24 in at 

10 GeV/c, which causes the required bend angle at to be 1.3 to 2.2 mrad greater 

than the nominal value. As a result, the actual momentum is expected to be l.TL 

below the value calculated from the nominal bend angle at at 20 GeV/c and 3.7% 

below this value at 10 GeV/c, when the quadrupoles are off. 

When the quadrupoles are turned on, the tuned current in was found 

to be higher, so it is clear that some steering in opposition to the bend of is 

being done by the quads. From the 20 GeV/c data in Table VA, it is seen that this 

2.5 

increase is about 20.2 “ 12% of the higher value when converted to momentum. 

Since the source is off-axis, the beam enters outside the axis and is 
bent towards the axis in the quads, which is in qualitative agreement with observa- 
tions. Using the position of the first collimator (0.64 in outside the axis) and 
the focal length of the Q5Q5 pair (200 in) we estimate that a ray could be bent 

through = 3.2 mrad. This is to be compared with the nominal bend angle 

3.2 

0, =59.3 mrad, giving = 5.47o. Thus after turning the quads on we might 

1 5y . 3 

expect to have to increase the current in by 5.47.. The fact that a larger 
increase is necessary suggests that our knowledge of the actual beam path prior 
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to remains ambiguous. This may be due to incomplete knowledge of the AGS 
fringing fields, or a misunderstanding as to the location of the off-axis 
components . 

Using the experimentally determined value of 127o for the quadrupole steer- 
ing and TL for the off-axis location of Ty we predict that the actual momentum 
is lAT, below the value calculated from the current in at 20 GeV/c when the 
quads are on. Since data on the amount of steering done by are not avail- 

able for lower momenta it is assumed that the same factor can be applied. 

F. Combined Momentum Calculation 

With the quads oif, at nominal momentum 20 GeV/c, the momentum as calculated 
from is 17.65 GeV/c minus 2.6%, which gives 17.2 GeV/c. The analysis based on 
the bend at D 2 gives 19.8 GeV/c minus (10.1 + 2.2)7o which equals (17.8 + 0.4)GeV/c. 
As our best value of the momentum, we quote the average, with error bars that 
encompass the full range of values within the uncertainties of either analysis 
alone. Thus the central value of beam momentum in this instance is (17.5 + 0.7) 
GeV/c, or 17.5 GeV/c with an uncertainty of + 4%,, 

With the quads on, at nominal momentum 20 GeV/c, the momentum as calculated 
from would then be 20,1 GeV/c minus 14%,- which gives 17.3 GeV/c. The analysis 
based on the bend at D 2 gives 19.8 GeV/c minus (10.1 + 2.2)7o which equals (17.8 
+ 0.4) GeV/c. Combining as before gives a central value of beam momentum of 
(17.55 + 0.65) GeV/c, or 17.6 GeV/c with an uncertainty of + 47.. This procedure 
is used in calculating the beam momentum in all experimental runs with the quads 
on. The results of the calculations are presented in Table Va. 

Table Vb then gives the central value of beam momentum for each data run. 

These values are uncertain by + 47,. The resolution of the momentum distribution 
is, coincidentally, 4% FWHM as discussed in Section VIlA above. 
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APPENDIX A 


Transmission of a Beam Telescope 

A beam telescope which has apertures of unequal wi^h is shown below; 



1 "2 

All rays which pass through aperture Aj^ at angles between the axis (zero degrees) 

and + * are transmitted, so the transmission function in this region is a con- 

A 

stant, T^. No rays pass through which travel at angles greater than iji with the 
axis, so that the transmission function drops to zero at + i|(g. Finally, at 
angles between and there is a linear fall-off in the transmission function. 
Altogether this leads to a trapezoidal transmission function for the telescope, 
as shown below. In the special case where the apertures are equal, “ 0 and 
T (ijf) reduces to a triangular shape. 



When a beam telescope is combined with a momentum-defining system, the trans- 
mission as a function of angle with respect to the beam axis is equal to the 
transmission as a function of the bend angle for rays constrained to pass through 
a point source. If the point source is replaced by a finite aperture, the trans- 
mission as a function of bend angle, T (0), is modified. The corners are rounded 
and the limiting angles of the "trapezoid" become larger as shown below. This 
effect is confined to a region within + 69^/2 of the corners, where 6®^ is the 
full angular width of the aperture in question. 


L 


Therefore, the transmission function of the second bend system in the test 
beam discussed here is determined as follows: 

(1) The scintillators and S 2 form a beam telescope with angles 


4-1 
2 (612) 


= 2.45 mrad 


4 + 1 

ijfp 4.07 mrad 

B 2 (612) 

(2) The entrance aperture of the second collimator C 2 approximates a point 
with angular width 

5 0 = = 1.07 mrad 

^ 950.5 - (460.75 - 24.0) 

(3) The shoulders of the transmission curve occur at 9=9- and 

9„ = 9„ + \jf.. Therefore, the telescope has a full width to the shoulders 


6 9^ = 2 i|r^ = 4,9 mrad 

This is to be compared to the nominal bend angle 62 ~ 24.2 mrad, and con- 
verted to momentum, thus: 


P2 . %-h - «H - Oj. . '>•’ . 


20X (Full Width to Shoulders) 


(4) The shoulders of the transmission curve are smeared out by + 

+ 0.54 mrad = + 2.2% due to the finite aperture of the collimator C 2 . 
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Preliminary Results 


The figures presented in this Appendix are produced 
from computer printout showing a preliminary analysis of 
typical calibration runs made with 20, 15, and 10 GeV 
protons. The events included were selected to have inter 
acted in one of the first 4 iron modules (2 interaction 
lengths). Figures B1 through B3 show the total number of 
particles, normalized to a mean of 100, detected in the 
first 7 iron modules (3-1/2 interaction lengths), for 20 
15, and 10 GeV incident protons respectively. 

Figures B4 through B6 show the number of particles, 
normalized to a mean of 100, detected at cascade maximum 
for 20, 15, and 10 GeV protons respectively. 
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TABLE 1 

layout of G-10 + 4.7° TEST BEAM, 18-21 July, 1970 


Element 


Symbol 

Location 

Description 

G-10 Target 


T 

Between G-10 & G-11 
ring magnets 

BeO wire, described in 
Reference (8) 

Virtual Target 



On a radial line 
through GIO target 

BeO wire, offset by fringe 
field of Gll 

Collimator #1 


^1 

268.25" from T 

1" hor X 4" vert x 36" lead 

Quadrupole #5 


Q5 

563" from T on 
"magnet axis" 

Type 8Q48 

Quadrupole #6 


^6 

623" from T on 
"magnet axis" 

Type 8Q48 

Bending Magnet 

#1 

°1 

701" from T 

Type 18D72 - bend angle 
3.4025° w.r.t. "magnet axis" 

Collimator #2 


^2 

460.75" from 

1/2" X 1/2" X 48" lead, with 
brass shims 

Collimator #3 


S 

535.75" from 

Variable width x 4" vert. 

X 48" lead automatic colli- 
mator 

Bending Magnet 

#2 

°2 

950.5" from 

type 18D72 - bend angle 
1.3916° 

Scint. Counter 

#1 


500" from I >2 
(Nominal) 

4" wide X 3-1/2" high plasti 

Cerenkov Counter 

C 

549" from 
(to entrance) 

2 meter long gas Cerenkov 
counter 

Scint. Counter 

#2 

^2 

612" from S^ 

1" X 1" plastic 

Scint. Counter 

#3 

^3 

55" from S2 

2" X 2" plastic 

NASA Experiment 

IS 

2258" from D 

2 

Ionization Spectrometer 



-26- 


TABLE II 
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TABLE IV 

Threshold Pressures for Cerenkov Counter 
Filled with Freon - 12 (CCI 2 F 2 ) 


(Calculated from n = l/P, and using Argonne National 

threshold 

Laboratory Report #6916, "Index and Dispersion of Cerenkov 
Counter Gases") 


PARTICLE 


MOMENTUM 
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(550) 
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17 
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20 GeV/c 

16 

5 

0.3 


All pressures in psia 
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(2) This table includes all data from situations when the beam magnets were actually tuned 
maximum flux. Many of these are not production runs. 



MOMENTUM BNL ANALOG DVM READINGS CORRECTED MOMENTUM 
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(1) The momentum calculation is described in the test. The best estimate of the central 
value is the combined value, where listed. This is known to about + 4%, while the 
resulution is, coincidentally, 4% FWHM. 
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Figure 2 


Figure 3 
Figure 4 


Schematic Cross Section of High Energy Cosmic Ray Experiment 

Schematic of G-10 + 4.7° Test-Beam Transport System, 

18 - 21 July, 1970 

High Energy Cosmic Ray Experiment (HECRE) Test Mounted 
Spark Chamber Orientation 
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